Spinor Bose-Einstein Condensates with Many Vortices 
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Vortex-lattice structures of antiferromagnetic spinor Bose-Einstein condensates with hyperfine 
spin F=l are investigated theoretically based on the Ginzburg-Pitaevskii equations near Tc. The 
Abrikosov lattice with clear core regions are found never stable at any rotation drive fi. Instead, 
each component ^l/i (j = 0,±l) prefers to shift the core locations from the others to realize almost 
uniform order-parameter amplitude with complicated magnetic-moment configurations. This system 
is characterized by many competing metastable structures so that quite a variety of vortices may 
be realized with a small change in external parameters. 
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Realizations of the Bose-Einstein condensation (BEC) 
in atomic gases have opened up a novel research field in 
quantized vortices as created recently with various tech- 
niques H, m ^ . Especially interesting in these systems 
are vortices of spinor BEC's in optically trapped ^^Na ||] 
and ^^Rb Q, where new physics absent in superconduc- 
tors §, "He g, and ^He §, |lO[ |ll|, may be found. 

Theoretical investigations on spinor BEC's were 
started by Ohmi and Machida and Ho followed 
by detailed studies on vortices with a single circulation 
quantum [l7[ 111, 111, |||, H, || . However, 

no calculations have been performed yet on structures in 
rapid rotation where the trap potential will play a less im- 
portant role. Indeed, the clear hexagonal-lattice image of 
magnetically trapped ^'^Na 1^ suggests that predictions 
on infinite systems are more appropriate for BEC's with 
many vortices. Such calculations have been carried out 
by Ho for the single-component BEC [p4[ and by Mueller 
and Ho for a two-component BEC [EST near the upper 
critical angular velocity ilc2 at T = 0. 

The purpose of the present paper is to perform de- 
tailed calculations on vortices of F = 1 spinor BEC's 
in rapid rotation to clarify their essential features. To 
this end, we focus on the mean-field high-density phase 
rather than the low-density correlated liquid phase [^ , 
and use the phenomenological Ginzburg-Pitaevskii (or 
Ginzburg-Landau) equations near Tc [p7| , p8| instead of 
the Gross-Pitaevskii equations at T = 0. Since fluctua- 
tions are small in the present system, this approach will 
yield quantitatively correct results near Tc- It should 
be noted that the corresponding free-energy is formally 
equivalent to that derived with Ho's "mean-field quan- 
tum Hall regime" near flc2 at T = ^ , so that the 
results obtained here are also applicable to that region. 

Model. — The free-energy density of an _F = 1 spinor 
BEC near Tc may be expanded with respect to the order 
parameters Vl/i {i = 0,±1) as 
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Here a, Pn , and (3s are expansion parameters, (fj, — 
X, y, z) denotes the spin operator, M is the particle mass, 
and summations over repeated indices are implied. The 
rotation axis is taken along z. The quantities /3n and ps 
are assumed to be constant near Tc with /?„ >0, whereas 
a changes its sign at Tc with a > for T < Tc. To sim- 
plify Eq. (^, we measure the length, the energy density, 
the angular velocity, and the order parameter in units 
of h/V2Ma, //3n, a/h, and a/(3n, respectively. The 
corresponding free-energy density is obtained from Eq. 
(Ill) by a ^ 1, /?„ 1, h^/2M 1, MVl and 
/3s — >■ <?s = /3s/ Pn- It thus takes a simple form with only two 
parameters (^s , ^) ■ We then introduce a couple of opera- 
tors by a = i{dx+idy) / \/2 and a'^ = £{—dx+idy)/\/2 with 
5= V — fJ^xr and which satisfy aa^— ata=l. 

Equation (m) now reads 



with rJc2 = l, and the free-energy is given by 



T= / /(r)dr. 
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We can find the stable structure for each {g^, fl) by min- 
imizing We have performed extensive calculations 
over the whole antiferromagnetic region gs > 0, where 
(^'i,*o,^'-i) = c'''Z^(0,l,0) and/ = i at 17 = with 6* an 
arbitrary phase and U the spin-space rotation | |l3| . 

A major difference from superfluid ^He lies in the 
fact that terms such as ^'*aa5'j [i^j] are absent, i.e., 
there are no gradient couplings between different compo- 
nents. Hence r2c2 is the same for all components, whereas 
in '^He only a single component becomes finite at f2c2 
to realize the polar state ||ll|. This degenerate feature 
is what characterizes the present system to bring many 
competing metastable structures, as seen below. 

Method. — We minimize Eq. (|3| ) with the Landau-level- 
expansion method (LLX) jl^, |29| by expanding the order 
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parameters as 
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*.(r) = Vv (4) 
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with N the Landau-level index, q the magnetic Bloch 
vector, and V the system volume. The basis functions 
{'4'Nq} are eigenstates of the magnetic translation oper- 
ator Tr = exp[— R • (V + X r)], which can describe 
periodic vortex structures of all kinds |2^. Its explicit 
expression is given by 

^jVq(r)= ^ ^i[qy(y+e^q^/2)+nai^iy+e^q^-naiy/2)/i^] 

n=-7Vf/2+l 

^^-ixy/2i^-{x~i^qy~nai^f /2l^ 

/ 2n£/a2 (x-f'qy-naij\ 
"V 2A-iV!0Fv H ^ 

with J\ff the number of the circulation quantum K = h/M 
in the system, Hj^ the Hermite polynomial, and the 
basic vectors in the xy plane with a2 || y and aixa2 = 
2t:P. Substituting Eq. (|) into Eq. (|) and carrying 
out the integration in terms of {s,t) = (x/aix,y/a2~ 
xaiy/aixa2), the free energy is transformed into a func- 
tional of the expansion coefficients {cjy'q}, the apex angle 
f3 = cos^^{aiy/ ax), and the ratio of the two basic vectors 
p=aila\ as = JF[{c^q}, /3, p]. For a given fi, we di- 
rectly minimize T jV with respect to these quantities. 

Search for stable structures. — We here sketch our 
strategy to find stable structures. To this end, we 
first summarize the basic features of the conventional 
Abrikosov lattice within the framework of LLX |^ : (i) 
any single q = qi suffices, due to the broken transla- 
tional symmetry of the vortex lattice; (ii) the triangular 
(square) lattice is made up of iV = 6n (4n) Landau levels 
(n = 0, 1, 2 • • • ); (iii) more general structures can be de- 
scribed by N — 2n levels, odd A^'s never mixing up since 
those bases have finite amplitudes at the cores. This 
Abrikosov lattice has a single circulation quantum per 
unit cell. 

With multi-component order-parameters, there can be 
a wide variety of vortices, which may be divided into 
two categories. We call the first category as "shift-core" 
states, where core locations are different among the three 
components with an enlarged unit cell. General struc- 
tures with Uk circulation quanta per unit cell can be de- 
scribed by using n,^ different q's, where the unit cell be- 
comes Uk times as large as that of the Abrikosov lattice. 
For example, structures with two quanta per unit cell are 
given by choosing (qi, q2) = (0, ^^~^^^ ), where bi and b2 
are reciprocal lattice vectors. This possibility was already 
considered by Mueller and Ho ||2^ for a two-component 
system and shown to yield stable structures. It also de- 
scribes the mixed-twist lattice to be found in •^He [ pT| . 



The second category may be called "fill-core" states with 
a single circulation quantum per unit cell (i.e., a single q 
is relevant). Here the cores of the conventional Abrikosov 
lattice are filled in by some superfluid components using 
odd-A^ wavefunctions of Eq. (||). This entry of odd-A^ 
Landau levels occurs as a second-order transition below 
some critical angular velocity smaller than flc2/S. It has 
been shown that the A-phase-core vortex in the _B-phase 
of superfluid ^He belong to this category 

We have carried out an extensive search for stable 
structures with up to ti^ = 9 circulation quanta per 
unit cell, including fill-core states. Since we are near 
Tc where normal particles are also present, we have per- 
formed the minimization without specifying the value 
of the magnetic moment M for the superfluid compo- 
nents. However, all the stable states found below have 
M = 0. Each of the three components were expanded 
as Eq. (H) using different q's, and the free-energy 
(^ is minimized with respect to c^^, /3, and p by us- 
ing Powell's method jso). To pick out stable structures 
correctly, we calculated Eq. ^) many times starting from 
different initial values for Cgq, /3, and p given randomly 

within 0.1 < RecJ,'^, Im4'^ < 0.9, O.Itt < f3 < O.Stt, and 
0.8 < /3 < 3, respectively. The state with the lowest en- 
ergy was thereby identified as the stable structure. The 
spin quantization axis and an overall phase were fixed 
conveniently to perform efficient calculations. Thus, any 
structures obtained from the solutions below with a spin- 
space rotation and a gauge transformation are also stable. 

Instability of Abrikosov lattice. — The present calcula- 
tions have revealed that the Abrikosov lattice with clear 
core regions is never stable at any rotation drive Q over 
the entire antiferromagnetic region gg > 0. Thus, any 
optical experiments to detect vortices by amplitude re- 
ductions are not suitable for the spinor vortices. 

The most stable Abrikosov lattice is given by 

*o(r) = Vy5]c^ViVq,(r), (6) 

JV 

with c^qj real, A^ = 6n, f-, and p=l. Here the an- 
tiferromagnetic component ^'g forms a hexagonal lattice 
with a single circulation quantum per unit cell. Below 
some critical velocity J7f smaller than Qc2/3, the core 
regions start to be filled in by 

^±,{r)^iVvJ2c^^UN'M^ (7) 

N' 

with A^' odd. The second transition for this odd-Landau- 
level entry occurs at ilf = 0.1497f2c2 and 0.0938ric2 for 
(7s = 0.1 and 0.3, respectively. 

However, calculations down to 0.0001ilc2 of using 1800 
Landau levels for gs>0 have clarified that the above fill- 
core state carries higher free energy than the following 



3 



0.202 



0.240 -0.02 



0.02 













1 


m 


o 


o 




1 








i 


1 


o 


o 


m 




1^ 












(a) 1^1 



(b)M^ 



FIG. 1: Spatial variations of (a) the order-parameter ampli- 
tude and rb) the magnetic moment along z, for the shift-core 
state (|)-(|) over |a;|, \y\ < ai at Qs = 0.08 and Q. = 0.95J7c2. 
The moment is directed along z. 



shift-core state with two circulation quanta per unit cell: 
^iir) = VvJ2cN^NciAr), (8) 

N 

4'_i(r) = Vy ^cw^wq,(r), (9) 

N 

with cn real and common to both, N = An, /3=-|, 
and (qi, q2) = (0, ^^^). The cores of *±i(r) are shifted 
from each other by ^(ai-|-a2). Figure 1 displays basic 
features of this shift-core state. The magnetic moment is 
ordered antiferromagnetically along z axis, and the am- 
plitude is almost uniform taking its maximum at each 
site where the moment vanishes to form the antiferro- 
magnetic state realized in the uniform state. 

Stable structures near nc2- — Having seen that the con- 
ventional Abrikosov lattice is never favored in the whole 
antiferromagnetic domain gs > 0, we now enumerate all 
the stable structures found near f2c2 to clarify their essen- 
tial features. This rapidly rotating domain is especially 
interesting, because the same free energy also becomes 
relevant near Qc2 at r = 0, as shown by Ho using a "mean- 
field quantum Hall regime" Thus, the conclusions 
obtained here are also applicable to the region at T = 0. 

Figure 2 displays the lowest free energy per unit volume 
as a function of gs for il. — 0.95ilc2- The value of each 
denotes the number of circulation quanta per unit cell. A 
special feature to be noted is that these various structures 
are energetically quite close to each other; for example, 
the = 8 state at gs = 0.02 is favored over the n„ = 3 
state by a relative free-energy difference of order 10^^. 
This fact suggests that we may realize quite a variety of 
metastable structures by a small change of the boundary 
conditions, the rotation process, etc. 

Details of these structures are summarized as follows. 

The = 2 state of Eqs. (^) and is stable for gs > 
0.0671. We have already seen its basic features above. 
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FIG. 2: Calculated free energy per unit volume as a function 
of gig at Sl = 0.95Slc2. Five different structures have been found 
for Qa > 0, and the value of each rin denotes the number of 
circulation quanta per unit cell for each stable structure. 



The n 1^ — 3,5,7 states can be expressed compactly as 
*o = E E W - e-'^^>^q„_„ ) , (10) 
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where A^'s are even and q^, — ;f-bi. These = 3,5,7 
states are stable for 0< 5s < 0.0196, 0.0313 < < 0.0613, 
and 0.0613 < 5s < 0.0671, respectively. Unlike the two 
component system considered by Mueller and Ho |2|] 
where each component is specified by a single-q basis 
function, 5'„i here is made up of multiple basis func- 
tions ■0A'q„ whose cores are shifted from each other by 
{i'/n^)a2- Figure 3 displays the basic features of the 

— 3 state. The lattice is hexagonal at 5s = 0, but de- 
forms rapidly as 5s increases. The order-parameter am- 
plitude is again almost constant, and the magnetic mo- 
ment M has a complicated structure. These features are 
common to all the > 3 states discussed here, although 
no details are presented for the other states. The lattice 
parameters {(}, p) for —5,7 states are (0.131l7r, 2.056) 
and (0.18307r, 2.609) for 5s = 0.05 and 0.066, respectively, 
which change little in each relevant range of stability. 

The remaining n„ = 8 state, stable over 0.0196 < 5s < 
0.0313, is given by 

*0 = X! '^Nqi (^Afqi + V'JVqs - "^Wqa + i'^Nc^) , (12) 
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FIG. 3: (a) Variations oi P/tt and p as a function of for the 
riK = 3 state at = 0.95r2c2 • This lattice for gia = is hexagonal 
with /3 = 7r/3 and p = 3. Figures 2(b)-(d) display, for (jis = 0.01, 
spatial variations of (b) the order-parameter amplitude, (c) 
amplitude of the magnetic moment, and (d) the magnetic 
moment, over — 3aii:/2 < a; < 3aii:/2 and —02/2 < y < 02/2. 
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0.025 are (O.SSITtt, 1.049), and 
changes only slightly in the above range of gs- 

Concluding remarks. — We have performed extensive 
calculations on antiferromagnetic F — I spinor vortices 
in rapid rotation. The conventional Abrikosov lattice is 
shown never favored. Each stable structure has almost 
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constant order-parameter amplitude and a complicated 



magnetic-moment configuration, as shown in Figs. 1 and 
3, for example. This means that any optical experiments 
to detect vortices by amplitude reduction will not be 
suitable for the spinor vortices. Instead, techniques to 
directly capture spatial magnetic-moment configurations 



will be required. The system has many mctastablc struc- 
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tures which are different in the number of circulation 
quanta per unit cell n„, but are quite close to each other 
energetically. Thus, domains to separate different struc- 
tures may be produced easily. This degenerate feature is 
also present within the (/3, p) space of a fixed n^, where 
is the vortex-lattice apex angle and p is the length ratio 
of the two basic vectors. Put it differently, we can de- 
form a stable lattice structure with a tiny cost of energy. 
These facts indicate that the spinor EEC's can be a rich 
source of novel vortices realized with a small change in 
external parameters. 
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